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Аннотация 

Введение. Безопасность судоходства и разработок подводных месторождений полезных ископаемых требуют 
точного обнаружения различных подводных объектов. В литературе рассматриваются вопросы отслеживания их 
перемещений и траектории движения. Предлагаются методы гидролокации, обеспечивающие высокую точность 
позиционирования подводных объектов. Отмечена высокая точность пеленга стереодатчиков с ультракороткой 
базой. Однако такое оборудование чувствительно к частоте дискретизации сигналов, что вызывает «шум дискре- 
тизации». В открытом доступе нет публикаций, посвященных решению этой проблемы. Представленное иссле- 
дование призвано восполнить данный пробел. Цель работы — изучение возможности получения данных, уточ- 
няющих информацию о пеленге подводных объектов за счет использования фазовой информации отраженных 
зондирующих сигналов и дополнительной процедуры передискретизации исходных данных. 

Материалы и методы. Местоположение объекта определяли с помощью экспериментального комплекса для 
исследования гидроакустических датчиков, созданного В.А. Широковым и В.Н. Милич в Удмуртском федераль- 
ном исследовательском центре Уральского отделения Российской академии наук. Использовали стереодатчик с 
малой базой (30 мм) по сравнению с расстоянием до объекта (=800-900 мм). Для обработки данных применяли 
методы цифровой фильтрации и математический аппарат корреляционного анализа отраженных гидроакустиче- 
ских сигналов, полученных фазовым методом. 

Результаты исследования. Представлены итоги сопоставления двух способов определения пеленга на объект: 
по разности времени прихода передних фронтов импульсов и по максимуму кросс-корреляционной функции 
(ККФ). Графически показано изменение пеленга при движении объекта. Использование переднего фронта сиг- 
нала обусловило небольшие выбросы значений вдоль всей кривой пеленга (менее 0,12 рад). При максимуме ККФ 
выбросы фиксировались лишь в некоторых областях, но были довольно значительными (около 0,17 рад). Пока- 
зано, как выбрать точки, соответствующие более гладкой и валидной траектории объекта, и как работать с оши- 
бочными точками. Представленный метод устранения ошибки можно реализовать программно. При квазигармо- 
ничном сигнале редкие измерения исходного сигнала интерполируются частыми вычисленными значениями. 
Благодаря такому виртуальному увеличению частоты дискретизации (передискретизации) можно фиксировать 
промежуточные показатели в оцифрованных исходных данных. Интерполяция значений сигнала кубическим 
сплайном позволила получить 20 точек на 1 период сигнала вместо 5 точек в исходном варианте. В этом случае 
более корректна траектория, сформированная с максимумом ККФ. 

Обсуждение и заключение. Задачу пеленгации можно решить с точностью, необходимой для практического 
применения. Учет фактора гладкости и непрерывности траектории движения объекта позволяет качественно кор- 
ректировать выбор максимума кросс-корреляционной функции сигналов стереодатчика. Предложенные методы 
обладают большим потенциалом для разработки систем подводного видения. 
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Abstract 

Introduction. Safety of navigation and development of underwater mineral deposits require the accurate detection of 
various underwater objects. The literature discusses the issues of tracking their motion and trajectory. Sonar methods are 
proposed to maintain high accuracy of underwater object positioning. High accuracy of the bearing of stereo sensors with 
an ultrashort base is noted. However, this equipment 1s sensitive to the sampling rate of the signals, which causes 
"sampling noise". There are no publicly available publications dedicated to the solution to this problem. The presented 
study is designed to fill this gap. This work is aimed to study the possibility of obtaining data clarifying information about 
the bearing of underwater objects through using the phase information about echoed probing signals and an additional 
procedure for resampling the source data. 

Materials and Methods. The location of the object was determined using the experimental complex for studying 
hydroacoustic sensors created by V.A. Shirokov and V.N. Milich at the Udmurt Federal Research Center, the Ural Branch 
of the Russian Academy of Sciences. A stereo sensor with a small base (30 mm) was used compared to the distance to 
the object («800—900 mm). Digital filtering methods and mathematical apparatus of correlation analysis of return 
hydroacoustic signals obtained by the phase method were used for data processing. 

Results. The results of comparing two methods for determining the bearing on an object are presented: by the difference 
in the time of arrival of the pulse-leading edges and by the maximum of the cross-correlation function (CCF). The change 
in bearing as the object moves, is graphically shown. The use of the leading edge of the signal caused small outliers of 
values along the entire bearing curve (less than 0.12 rad). At the maximum CCF, emissions were recorded only in some 
areas, but they were quite significant (about 0.17 rad). It showed how to select points corresponding to a smoother and 
more valid object trajectory, and how to work with erroneous points. The presented method of error correction can be 
implemented programmatically. With a quasi-harmonious signal, rare measurements of the original signal are interpolated 
by frequent calculated values. Thanks to this virtual increase in the sampling rate (oversampling), intermediate indicators 
can be recorded in the digitized source data. Interpolation of the signal values by a cubic spline allowed us to obtain 20 
points for 1 period of the signal instead of 5 points in the original version. In this case, the trajectory formed with the 
maximum CCF is more correct. 

Discussion and Conclusion. The direction-finding problem can be solved with the accuracy required for practical 
application. Taking into account the factor of smoothness and continuity of the object's trajectory makes it possible to 
qualitatively correct the selection of the maximum of the cross-correlation function of the stereo sensor signals. The 
proposed methods have great potential for the development of underwater vision systems. 


Keywords: determination of the location of an object in a hydraulic environment, phase direction finding, pulse-leading 
edges, cross-correlation function, sampling noise 
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Введение. Обеспечение безопасности судоходства и работы на подводных месторождениях полезных HCKO- 
паемых требуют качественного обнаружения подводных объектов и отслеживания их перемещений [1]. В систе- 
мах подводного наблюдения применяются гидроакустические датчики. Они улавливают сигнал, отраженный от 
объекта, и позволяют вычислить его местоположение методом трилатерации [2]. В этом процессе каждый датчик 
предоставляет информацию об интервале времени прохождения зондирующего сигнала, отраженного от объекта. 
В случае использования нескольких датчиков [3] возникает возможность решения задачи обратной простран- 
ственной засечки и определения координат наблюдаемого объекта [4]. Для повышения точности измерений пред- 
ставляется перспективным использование в качестве приемника стереодатчиков с ультракороткой базой [5], поз- 
воляющих получать фазовую информацию [6] и определять пеленг на объект [7]. Изучены возможности исполь- 
зовании гидролокационных данных об удаленных подводных целях, а также о пеленге подводных объектов [8] 
за счет использования фазовой [9] или частотной [10] информации отраженных зондирующих сигналов. 

При этом известно, что процедура дискретизации сигнала обусловливает ошибки при определении парамет- 
ров траектории движения подводных объектов [11]. Это относится как к дальномерным измерениям, так и к про- 
цедурам определения пеленга. Оцифровка аналогового сигнала датчика, необходимая для дальнейшей цифровой 
обработки, вносит так называемую «погрешность дискретизации». Погрешность, вызванная квантованием ам- 
плитуды сигнала, оценивается количеством уровней квантования аналого-цифрового преобразователя. Погреш- 
ность, вызванная дискретизацией по времени, пропорциональна величине интервала времени квантования и ско- 
рости изменения напряжения. Поэтому частоту дискретизации сигнала стремятся сделать как можно выше верх- 
ней частоты самого сигнала. Однако увеличение частоты дискретизации во многих случаях ограничивается воз- 
можностями регистрирующей аппаратуры: процессора, аналого-цифрового преобразователя, каналов передачи 
данных, устройствами хранения данных. Ограниченность возможностей дискретизации во времени при выпол- 
нении траекторных измерений не позволяет точно фиксировать экстремальные значения траектории (дальность 
и пеленг). Как следствие, снижается точность результатов измерений, появляются выбросы на фиксируемых тра- 
екториях. Для устранения этого недостатка необходимо изучить потенциал передискретизации оцифрованного 
гидроакустического сигнала. В открытом доступе нет публикаций, посвященных решению данной проблемы. 
Имеющийся пробел восполняют материалы этой статьи. 

Цель представленного исследования состоит в оценке возможностей определения пеленга на объект с исполь- 
зованием фазовой информации и в отработке фазового метода пеленгования с применением передискретизации 
оцифрованного гидроакустического сигнала. 

Материалы и методы. Рассматривается алгоритм определения координат с применением пеленгации фазо- 
вым методом. Приводятся результаты его тестирования для объекта, движущегося по круговой траектории в экс- 
периментальном бассейне. 

Постановка эксперимента. При проведении опытов для определения местоположения объекта в гидросреде 
при пеленгации фазовым методом задействовали экспериментальную установку, созданную В.А. Широковым и 
В.Н. Милич в Удмуртском федеральном исследовательском центре Уральского отделения Российской академии 
наук [12]. Это лабораторный измерительный комплекс с линейной аквасредой в виде протяженного цилиндриче- 
ского резервуара (гидроволновода) и бассейном, оснащенным системой генерации испытательных гидроакусти- 
ческих сигналов. Ниже детально описаны элементы установки, которые использовались в опытах. 

1. Излучатель гидроакустического сигнала с известными координатами 5 (рис. 1). В работе используется ам- 
плитудно-модулированный сигнал [13]. 


M 


Рис. 1. Расположение принимающих датчиков (стереодатчик с двумя приемниками 51 и 512) и излучателя (5) в плоскости 
экспериментального бассейна (xy). Здесь d — база стереодатчика; D(O; st) — расстояние между стереодатчиком и объектом; 
а — пеленг на объект; плоскость (x у’) образована плоскостью стереодатчика и нормалью и к ней 
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2. Отражающий объект. В эксперименте в качестве тестового используется протяженный цилиндрический 
объект — медный провод, диаметр которого (0,15 мм) значительно меньше длины акустической волны (1,5 мм). 
В представленной работе задача решается на плоскости, поэтому для составления вычислительных процедур 
расчета координат данный объект можно считать точечным. 

3. Стереодатчик, преобразующий гидроакустический сигнал в электрический. Известны координаты каждого 
из двух приемников стереодатчика (511, 512 на рис. 1). 

4. Аппаратно-программный комплекс, выполняющий усиление, оцифровку и обработку сигналов датчиков. 

Результатом измерений является файл регистрации принимаемого двухканального сигнала. 

Алгоритм обработки результатов измерений. Использование стереодатчика с малой базой (30 мм) по срав- 
нению с расстоянием до объекта (= 800—900 мм) позволяет применить упрощенные формулы. 

Итак, рассчитаем координаты исследуемого объекта в системе x'y', образованной плоскостью стереодатчика 
и нормалью к ней. Для этого задействуем алгоритм, включающий пять этапов. 

1. Предварительная обработка: 

— удаление слепой зоны (первые М отсчетов); 

— фильтрация сигнала. 

В работе применяется высокочастотная фильтрация. 

2. Выделение из всего сигнала информативных фрагментов impl, imp2 (в двух каналах), содержащих им- 
пульсы, отраженные от исследуемого объекта. Для этого выполняется пороговая обработка по амплитуде А сиг- 
нала. Порог принимается равным рхтах(А) (р = 0,5). 

Такая обработка позволяет определить передний фронт импульса (BeginIndex), отраженного от объекта. Из двух 
каналов сигнала вырезаются участки [BeginIndex-0,5xLenSignal, BeginIndex*LenSignal]. Здесь LenSignal — длина car- 
нала, подаваемого на вход. В результате получаем два сигнала impl и ітр2 длиной 1,5xLenSignal. В каждом из них 
содержится импульс, отраженный от объекта. 

3. Расчет расстояния D(O; sf) между объектом и стереодатчиком. 

— Определение расстояния, пройденного импульсом от излучателя до объекта и от объекта до датчика при- 
ема, по формуле D(S,O,st)-BeginIndexxdtxC. Здесь С=1475 м/с — скорость распространения акустической волны, 
dt = 0,2x10—6 c — интервал дисктретизации. 

— Расчет значения: 

D( S,O,st )— D( S,st ) 


D(O,st)- 3 (1) 
4. Определение пеленга а на объект [14]: 
a = arcsin Ат (2) 


где À — длина волны; F — частота сигнала; dt — интервал дискретизации; т — разница прихода отраженного 
импульса на два стереодатчика (в отсчетах); d — база стереодатчика (в представленном исследовании À = 1,5 MM, 
Е= 1 МГЦ, dt = 0,2 мкс, d = 30 мм). 

5. Вычисление координат анализируемого объекта в системе, образованной плоскостью стереодатчика и его 
нормалью (рис. 1), по формулам: x —D(O, st) *cos(a); y =D(O, sf) xsin(a). 

Определение пеленга на объект. Примем условие малости стереобазы по сравнению с расстоянием OT дат- 
чика до объекта. В этом случае стереодатчик позволяет определить пеленг на объект с помощью информации о 
разности фаз принятых импульсов, ранее отраженных от объекта, а также по разности времени прихода передних 
фронтов импульсов, поступивших на стереодатчик (рис. 2). 


A, B, 


Th1 
Th2 


Рис. 2. Определение разности моментов прихода ds стереосигналов, которые обозначены сплошной imp 1 (f) 
и пунктирной imp2(t) линиями. Thl, Th2 — пороги для первого и второго сигнала стереопары соответственно; 
dt — интервал дискретизации; А — значение сигнала B вольтах; f — время 
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Для реализации такого подхода выполняются описанные ниже процедуры. 

1. Идентификация моментов прихода передних фронтов стереосигналов с помощью пороговой обработки [15] 
с порогом рхтах(А). В данной работе коэффициент р принят равным 0,5 [16]. 

2. Вычисление разности ds (количество интервалов df) между передними фронтами импульсов, зарегистриро- 
ванных в двух приемниках стереодатчика. 

3. Вычисление пеленга а по формуле (1). За т принимаем значение ds. 

При таком способе вычисления сдвига фазы передний фронт сигнала будет главным индикатором импульса, 
отраженного от объекта. Не учитывается внутренняя структура сигнала, которая может иметь решающее значе- 
ние для корректного определения пеленга на объект. Учтем схожесть принимаемых сигналов и рассмотрим в 
работе другой подход к вычислению пеленга — на основе кросс-корреляционной функции (ККФ) стереосигна- 
лов. Логично определять сдвиг фазы по максимуму ККФ. Ниже описаны элементы алгоритма вычисления пе- 
ленга с помощью ККФ. 

1. Расчет ККФ импульсов, поступивших на два приемника одного стереодатчика: г = xcorr(impl, imp2). Зна- 
чение функции xcorr при смещении m второго сигнала относительно первого вычисляется следующим образом: 
N-m-l 

impl,,,,x imp2,, m 2 0, 
XCOIT;,51 imp2 (т) = =, 
ХСОГТ2 тр (7), m < 0. 

Здесь № — длина impl и imp2, -N«m«N. 

2. Определение максимума ККФ и соответствующего ему сдвига dc (в orcuerax). 

3. Вычисление пеленга а по формуле (1) с учетом, что за т принимается величина dc. 

Результаты исследования 

Тестирование методов определения пеленга. Представленные способы определения пеленга на объект (по 
разности времени прихода передних фронтов импульсов и по максимуму ККФ) протестировали для объекта, 
движущегося по круговой траектории в экспериментальном бассейне. Использовался гармонический сигнал дли- 
тельностью 7 периодов и частотой 1 МГц. Оцифрованные стереосигналы регистрировались в 256 точках траек- 
тории движения объекта. Частота оцифровки — 5 МГц. 

На рис. 3 показано изменение значений пеленга по мере движения объекта по окружности для двух предло- 
женных подходов (разница передних фронтов, максимум ККФ). 
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Рис. 3. Изменение пеленга а на объект, совершающий круговое движение. Показатели получены при использовании двух 
подходов (красный маркер — разница времени прихода передних фронтов стереосигналов, синий маркер — максимум 
ККФ). MN — номер измерения, а — пеленг на объект 
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При использовании переднего фронта сигнала (рис. 3, красный маркер) почти вдоль всей кривой изменения 
пеленга на объект наблюдаются небольшие выбросы значений. При использовании максимума ККФ (рис. 3, си- 
ний маркер) только в некоторых областях присутствуют выбросы значений пеленга, но они достаточно большие. 
Возможно, это объясняется особенностями движения объекта в области, где зондирующий и отраженный сиг- 
налы пересекают зону турбулентности, вызванной движением объекта. При нахождении фазового сдвига стерео- 
сигналов с помощью максимума ККФ возникают дополнительные ошибки, что затем отражается в вычисленной 
траектории объекта. Также выбросы в вычисленных значениях пеленга могут быть обусловлены недостаточно 
высокой частотой дискретизации. Значения этих ошибок существенно выше показателей выбросов, наблюдае- 
мых для метода с использованием переднего фронта импульса. 

На рис. 4 а представлены траектории движения объекта, вычисленные с использованием глобального макси- 
мума ККФ и его соседних локальных максимумов. 
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Рис. 4. Анализ траекторий кругового движения объекта: а — траектории, вычисленные с использованием глобального 
максимума ККФ (черный маркер), локального ближайшего правого максимума ККФ (оранжевый маркер), локального 
ближайшего левого максимума ККФ (зеленый маркер); 6 — пример ККФ, глобальный максимум которой соответствует 
корректному значению координат объекта; в — пример ККФ, в которой корректному значению координат объекта 
соответствует левый локальный максимум ККФ 


Используя априорную информацию о плавности перемещения экспериментального объекта, можно выбрать 
точки, которые соответствуют более гладкой и валидной траектории. Такой подход применяется в других рабо- 
тах [17]. На рис. 4 а при х '€[800;870], y ^0 корректные значения находятся на кривой для траектории, рассчитан- 
ной по левому локальному максимуму ККФ. Есть также ошибочные точки. Если их заменить точками траекто- 
рии, полученной в результате использования правого локального максимума ККФ, можно вычислить корректное 
значение точки траектории (x 800, у’<0 на рис. 4 а). Данный графический метод устранения ошибки может быть 
реализован программно. 

Таким образом, в некоторых точках траектории корректный сдвиг фазы стереосигналов, необходимый для 
вычисления пеленга, может соответствовать не глобальному максимуму ККФ, а одному из соседних локальных 
экстремумов (рис. 4 6, в). Такое смещение максимумов ККФ может быть обусловлено недостаточно высокой ча- 
стотой дискретизации сигналов. Этим же объясняется ступенчатый характер результирующей траектории. 

С учетом квазигармоничности сигнала предлагается использовать интерполяцию редких измерений исход- 
ного сигнала частыми вычисленными значениями. Такое виртуальное увеличение частоты дискретизации 
(передискретизации) дает возможность зафиксировать промежуточные значения в оцифрованных исходных дан- 
ных. Интерполяция значений сигнала кубическим сплайном позволила получить 20 точек на 1 период сигнала 
вместо 5 точек в исходном варианте. На рис. 5 приведены расчетные траектории движения объекта, сформиро- 
ванные с предварительной передискретизацией сигналов. 
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Рис. 5. Расчетные траектории движения объекта, полученные для стереосигналов с виртуально повышенной частотой 
дискретизации: синий маркер — используется разница времени прихода передних фронтов сигналов; 
красный маркер — используется максимум ККФ) 


—30 1 


Красная линия — это траектория движения объекта, полученная c использованием максимума ККФ для pac- 
чета пеленга по предварительно интерполированным данным. Она представляется более корректной по сравне- 
нию с аналогичной траекторией, полученной с использованием разницы передних фронтов отраженных импуль- 
сов (синяя линия). 

Обсуждение и заключение. Анализ особенностей определения пеленга на подводный объект с использова- 
нием взаимной фазовой информации сигналов дифференциального стереодатчика дает возможность сделать сле- 
дующие выводы: 

— разность времени прихода передних фронтов импульсов, зарегистрированных стереодатчиком, позволяет получать 
значения пеленга, показывающие направление на объект. Качество результата зависит от принятой величины порога (ис- 
пользуемого при определении переднего фронта импульса) и изменчивости амплитуд принимаемых сигналов; 

— применение подхода с использованием максимума кросс-корреляционной функции при недостаточной ча- 
стоте дискретизации исходных сигналов приводит к существенным выбросам в расчетной траектории (рис. 3); 

— передискретизация сигналов с помощью интерполяции исходных квазигармонических сигналов с исполь- 
зованием максимума ККФ для расчета пеленга и координат объекта позволяет получить гладкую траекторию 
(рис. 5), соответствующую плавному перемещению объекта по кругу. Такой подход требует значительно боль- 
ших вычислительных затрат в сравнении с другими, рассмотренными в этой статье. Это может повлиять на ско- 
рость обработки сигналов в режиме реального времени. 

Таким образом, исследованные возможности уточнения фазового метода пеленгации (передискретизация 
оцифрованных сигналов и использование максимума ККФ) позволяют эффективно оценивать траекторию дви- 
жения подводного объекта. 

Описаны условия получения гладкой корректной траектории движения отслеживаемого объекта. Этой цели 
служит передискретизация оцифрованных с недостаточной частотой дискретизации интерполированных квази- 
гармонических стереосигналов. Вместе с этим задействуется метод расчета пеленга, использующий максимум 
кросс-корреляционной функции сигналов стереодатчика. 

Данные на рис. 5 подтверждают, что задачу пеленгации можно решить с точностью, необходимой для прак- 
тического применения. Учет фактора гладкости и непрерывности траектории движения объекта позволяет одно- 
значно корректировать выбор максимума кросс-корреляционной функции сигналов стереодатчика (рис. 4). 

Предложенные методы имеют большое значение для разработки систем подводного видения. 
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